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Abstract 

Battery electrode surfaces are generally coated with electronically insulating solid films of thick¬ 
ness 1-50 nm. Both electrons and LK can move at the electrode-surface film interface in response 
to the voltage, which adds complexity to the “electric double layer” (EDL). We apply Density 
Functional Theory (DFT) to investigate how the applied voltage is manifested as changes in the 
EDL at atomic lengthscales, including charge separation and interfacial dipole moments. Illus¬ 
trating examples include Li3P04, Li2C03, and Li3;Mn204 thin-films on Au(lll) surfaces under 
ultrahigh vacuum conditions. Adsorbed organic solvent molecules can strongly reduce voltages 
predicted in vacuum. We propose that manipulating surface dipoles, seldom discussed in battery 
studies, may be a viable strategy to improve electrode passivation. We also distinguish the com¬ 
puted potential governing electrons, which is the actual or instantaneous voltage, and the “lithium 
cohesive energy”-based voltage governing Li content widely reported in DFT calculations, which 
is a slower-responding self-consistency criterion at interfaces. This distinction is critical for a com¬ 
prehensive description of electrochemical activities on electrode surfaces, including Li"*" insertion 
dynamics, parasitic electrolyte decomposition, and electrodeposition at overpotentials. 


keywords: lithium ion batteries; voltage prediction; density functional theory; computational elec¬ 
trochemistry 


1 


I. INTRODUCTION 


Unlike pristine noble metal or graphite basal-plane electrodes used in classical electric 
double layer (EDL) studies,lithium ion battery (LIB) electrodes generally exhibit complex 
interfaces.-'^ Both electron (e“) and Li’*' transport can occur inside LIB electrodes. In 
addition, solid thin films, on the order 1-50 nm thick, are ubiquitous on LIB electrode 
surfaces and can become part of the EDL (Fig.[T^-h). The EDL is critical because key battery 
processes like Li’*' incorporation kinetics (including Li"*" desolvation), parasitic reactions, Li- 
plating, and degradation-inducing phase transformations on electrode surfaces most likely 
initiate within it. Indeed, “solid electrolyte interphase” (SEI) films^i^^- covering anodes 
are relied upon as kinetic barriers that prevent thermodynamic equilibrium in batteries; 
they are critical for LIB which operate outside the redox stability window of the organic 
solvent electrolytes used. As discussed below, the EDL is also crucial in computational work 
because it effectively determines the electrode potential (henceforth “applied voltage” or 
simply “voltage”) in the simulation cell. Our work examines EDL associated with thin-hlm 
coated electrode surfaces, and emphasizes the importance of creating electrode interface 
models that exhibit consistent electrochemical activities for both Li+ and e~. 

Examples of solid films on electrode surfaces include Li 2 C 03 layers formed on pris¬ 
tine cathode oxide surfaces;-'— cathode-coating Elms made of electrolyte decomposition 
products;-'^'— SEI films on anodes arising from reductive decomposition of liquid electrolyte 
components;-'^'- artificial protective/passivating coatings,— including atomic layer deposi¬ 
tion (ALD) layers^'^^ which can undergo phase transformations at low voltages;— ALD 
layers between solid electrolytes and electrodes in all-solid state batteries;— and even Li 202 
hlms deposited on cathodes during fast discharge of Li-air batteries,—'^ the re-oxidation of 
which is accompanied by signihcant voltage hysteresis and is a root cause of the lack of a 
suitable liquid electrolyte in Li-air batteries. 

Detailed atomic lengthscale understanding of the interfaces and EDL associated with 
such thin films has remained elusive, partly due to challenges in imaging buried interfaces 
at sub-nanometer resolution. Even EDLs associated with liquid organic electrolyte in bat¬ 
teries have arguably only recently received systematic experimental studies.—'— Modeling 
efforts have made much progress in elucidating the structures of solid-solid^^i'— and solid- 
vacuum^'— interfaces. However, voltage dependences, particularly in processes involving e~ 
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FIG. 1: Some scenarios for electric double layers (EDL). (a) Voltage drops at electrode/liquid 
electrolyte interface when the metallic electrode is uncharged. Organic solvent molecules orient 
themselves to reduce the vacuum voltage (Sec. lIVPp . (b) Solid films with fixed dipoles, like fer- 
roelectrics, can serve the same purpose as molecules, and represent a novel strategy for electrode 
passivation fSec. lIVBll . (c)-(h) are additions to intrinsic dipole-induced voltage drops of (a) when 
charges exist on the cathode active material and/or thin films, (c) Pristine, noble metal electrode 
in contact with liquid electrolyte, (d) Inert thin film intervening between metal and liquid, (e) 
EDL entirely at metallic electrode/thin film interface fSec. IIV "Ai IIV Cp . (f) Same as (e), with more 
charge separation (Sec. IIV Al) . (g) EDL is only partially inside thin film (Sec. IIVDI) . (h) EDL is 
between redox-active film and liquid electrolyte fSec. lIV^ . Yellow and red backgrounds depict Au 
and liquid electrolyte; light brown textures are solid electrolyte thin hlms or redox-active cathode 
materials; green circles are negatively charged counter ions in the liquid electrolyte. For simplicity, 
cations in the liquid are omitted. Voltage profiles and surface charge densities are caricatures; in 
the calculations, voltage changes are 1-3 V and film thicknesses (demarcated by vertical dashed 
lines) are about 1 nm. The plots represent static conditions; ohmic losses add slopes to all flat 
portions. 
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transfer like parasitic reactions and nndesirable metal plating, have arguably received less 
theoretical attention. In this work, we apply electronic Density Functional Theory (DFT) 
calculations to investigate how voltages affect the structures and stoichiometries of thin solid 
films, their interfaces, and their EDL at atomic lengthscales. In view of the complexities of 
LIB electrodes, we have adopted simple model systems. The liquid electrolyte is omitted, 
although a few solvent molecules are included as frozen monolayers in some models to illus¬ 
trate their huge impact on the voltage. Au(lll) surfaces are adopted as model cathodes, 
instead of transition metal oxides typically found in batteries. An does not alloy with Li un¬ 
der voltages considered in this work and is a convenient inert electrode. These systems might 
be realized in ultra-high vacuum (UHV) settings; they dovetail with the use of copper— and 
gold^ electrodes in recent fundamental battery science experimental studies. 

The model thin Elms examined in this work, in order of increasing complexity, are Li 3 P 04 
(010) (Fig.[2^),— Li 2 C 03 basal plane (Fig.[2)D),— and Lij.Mn 204 (111) (Fig.[2l3).— These 
are coated on Au(lll) on one side and face a vacuum region on the other. Li 3 P 04 is a non- 
redox-active solid electrolyte. It illustrates the correlation between interfacial dipole densi¬ 
ties and voltages. Li 2 C 03 is often found on as-synthesized cathode surfaces. In LIB studies, 
there is disagreement^ concerning whether Li 2 C 03 dissolves upon soaking in electro lyte,-^^ 
is formed at high voltages,— or is removed at >4 V.- In Li-air batteries, Li 2 C 03 is gener¬ 
ally accepted to be oxidized and removed above ~4.5 V.— Our goal is not to elucidate the 
detailed reaction mechanism, but to study the electronic and structural signatures accompa¬ 
nying voltage increase. Redox-active LiMn 204 (spinel LMO) is not used as nanometer-thick 
hlms in batteries, but its inclusion here helps illustrate the diversity of interfacial behavior 
(Fig. [1]). Like most LIB cathode materials, LMO is a polaron conductor because Mn can 
exhibit multiple redox states. It also conducts Li'*’, but is not a band (“metallic”) conduc¬ 
tor. By classical electrostatics dehnition, the interior of LMO is not instantaneously at a 
constant potential.— Our LMO/Au interface model provides a well-defined Fermi level {Ep) 
to demonstrate how voltage control is achieved at atomic lengthscales, and how this affects 
EC oxidative decomposition on LMO surfaces.— 

Out of necessity, we critically examine the voltage calibration method used in the theo- 
rtical LIB literature (Sec. [TTl)- The same experimental voltage governs Li'*' and e~ motion. 
Indeed, in cyclic voltametry, it is difficult to distinguish currents arising from Faradaic Li- 
insertion or parasitic processes that involve undesired e~ transfer to the liquid. But what 
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(a) (b) (c) 


FIG. 2: (a) Li 3 P 04 (010)/Au(lll). (b) Li 2 C 03 basal plane/Au(lll). (c) 

Lii_xMn 204 (lll)/Au(lll). Li, P, O, C, H, Mn, and An atoms are depicted in blue, lime- 
green, red, grey, white, purple, and yellow, respectively. 

is “voltage” at atomic lengthscales? The absolute electrochemical potential for a charged 
particle in an atomic conhguration is the (free) energy needed to transfer it into a given 
phase from inhnity (see, e.g.. Ref. ? ). There can be differences whether the particle is a 
Li’*' or a e“. In DFT modeling, the electronic voltage (Vg) governing electronic motion, and 
ionic voltage (Vg) governing Li'*' motion, can and must be distinguished. 

First we focus on Vg. Under constant voltage experimental conditions, potentiostats 
directly control the e~ (not Li) content. When the phase is a metallic conductor with a 
well-dehned Up, the voltage asociated with transferring e~ from inhnity, multiplied by the 
electronic charge, is equal to the work function (<F) modulo a constant, just like in photo- 
voltaics (PV)^ and electrocatalytic systems.— <F=Uvac — where is the vacuum level 
and E-p is set to zero herein without loss of generality. If the electrode is immersed in liquid 
electrolyte, $ calculated for the electrode covered by a sufficiently thick liquid layer, divied 
by the electronic charge, also gives Vg.— Vg applies even under out-of-equilibrium condi¬ 
tions. It is readily computed in our model systems because vacuum regions exist and qualify 
as “inhnity.” In commercial batteries, there is no vacuum. Nevertheless, our focus on <F is 
formally correct and emphasizes the synergy between batteries, fuel cells,— PV,— and 
metal-semiconductor interfaces.—*^ In those research areas, dipoles at interfaces are recog¬ 
nized as properties of paramount importance. At battery liquid-solid interfaces, EDL’s are 
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FIG. 3: Overview of representive predictions of how EDL affects electronic voltage (Ve). The white 
background is vacuum. Yellow and textured regions are Au and solid thin films, respectively: (a) 
Li 3 P 04 /Au fSec. UVAD : (b) LiaCOg/Au ISec. HVTTIl : (c) EC+PFg/LiaCOs/Au ISec. HVDll : (d) 
EC/Li 2 ;Mna 04 /Au (Sec. lIVEp . are negatively charged Li+ vacancies. Ve values with and 

without parentheses are two alternate scenarios. Underlined values are at equilibrium (i.e., Ve=Vi). 
PFg anions on Li 3 ;Mna 04 surfaces are found to accelerate EC decomposition in their vicinity (not 
shown; Sec. IIVEI) . 

the manifestations of complex surface dipole distributions associated with applied voltages; 
it is in effect our computational potentiostat. Based on our calculations, we propose the 
manipulation of surface dipoles as a novel electrode passivating strategy (Fig. (Hd). 

Li'*’ ions are only indirectly controlled by potentiostats; they slowly redistribute in re¬ 
sponse to changes in electronic configurations. We define the Li “voltage” (Vj) as the (free) 
energy difference between Li chemical potential in a system and the lithium metal cohesive 
energy per Li divided by |e|. This definition is widely used as the sole voltage estimate 
in the LIB theory literature.—'"— At equilibrium, when the voltage is pinned at the redox 
potential of Li-insertion reactions, Ve must be equal to Vj. Thus Vi should be considered a 
self-consistent criterion at interfaces. By itself, Vi can give incorrect voltages when model¬ 
ing parasitic processes induced by e~ transfer, and it may not govern the dynamics of Li'*' 
insertion processes adequately. 
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Fig. |3^-d offers a preview of our results, pertinent to issues discussed above. The rest of 
the paper is organized as follows. Sec. 2 discusses and reconciles different dehnitions of volt¬ 
age in DFT calculations in more detail. Sec. 3 describes the methods used. The results are 
given in Sec. 4 and Sec. 5 summarizes the paper. A supporting information (S.I.) document 
provides further details about our three model systems, compares the modeling of battery 
interfaces with geochemical^^— and other liquid-oxide interfaces,—!^ and discusses Li 2 C 03 
oxidation thermodynamics. Experimentalist readers are encouraged to skip to Sec. IIVI 


II. VOLTAGE CALIBRATION IN DFT CALCULATIONS 

DFT calculations are conducted at constant number of electrons, not constant electro¬ 
chemical potential. Estimating and maintaining voltages in periodically replicated, con¬ 
densed phase DFT simulation cells have long been recognized as a challenge.— Re¬ 
cently we have applied an ab initio molecular dynamics (AIMD)-based Li’*' transfer free 
energy method to predict voltage dependences of redox processes at interfaces between liq¬ 
uid electrolytes and lithium-intercalated graphite.—!^ AIMD has also been coupled with 
continuum approximations of liquid electrolytes to model constant voltage conditions on 
pristine electrodes.— But present-day AIMD simulation time scales are too short to extend 
to Li’*' motion inside LIB electrodes or the surface Elms coating them except at significantly 
elevated temperatures. This work focuses on how applied voltages affect solid Elms. 

Thus we model solid interfaces at zero temperature. With sufEcient equilibration time and 
in the presence of Li'*' and e~ reservoirs, the total number of electrons (ue) and Li’*" content 
(nLi+) in an electrode in a LIB are governed by the chemical potentials for electrons {fig) and 
Li"*" ions (/iLi+), respectively. E-p of a metallic electrode is /ie modulo a constant. Here we 
have not used fie or fiu+, the notations for total electrochemical potentials of e~ and Li"'“,T 
precisely because constant-charge DFT calculations and terminologies are being applied to 
infer constant-potential electrochemical properties. Surface potential contributions to /i’s 
should naturally be present in our simulation cells, which feature explicit interfaces. 


7 



A. Electronic Voltage (Ve) from Work Function 


The electronic voltage of an atomic configuration which is band conductor, imposed by a 
potentiostat, not necessarily at equilibrium with respect to nu or the Li+ occupation sites, 
is defined as 

Ve = $/|e|-1.37 V, (1) 


if the Li”^/Li(s) reference is used.— In the literature, values of 1.39 V and 1.44 V have 
been adopted; the small discrepancy is a measure of possible systematic error. In the DFT 
formulation, the e~ conhguration hnds its ground state instantaneously. Ve pertains to 


inhnitesimal changes in Ue, which leaves the system uncharged. Note that Refs. l72 and 
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deal with Ve despite the fact that Li^ is transferred, 
metallic electrode at its Fermi level. As stated in Ref. 


Decause the excess e is left on the 


73l . the free energy of that graphite 


edge-plane system has yet to be optimized with respect to the surface Li-content. In the 
present work, Ve are always reported for locally optimized atomic conhgurations so that all 
forces are zero. 


B. Lithium Metal Cohesive Energy (Vj) for Voltage 

In the “lithium cohesive energy” (Vj) approach, discrete and matched numbers of e~ and 
Li+ are added/subtracted simultaneously: 

PLi = E{nu+, rie) - E{nu+ - 1, Ue - 1) 

= E{nu+, Ue) - E{nu+ - 1, Ue) E{nu+, n^) - E{nu+: - 1) 

~ h'Li-l- T /^e- (2) 

The last equality holds when the system approaches inhnite size. Here E{nu+,ne) is the 
zero temperature total energy after optimization of all interior degrees of freedom. The 
T = 0 K voltage is then given by 


Vi = (/XLi - ^Li(s))/|e|, (3) 

where F^Li(s) is the cohesive energy of lithium metal, predicted to be 1.57 eV using the 
DFT/PBE functional. No net charge is introduced. There is an ambiguity concerning 
whether fiu{nu) should be E{nu -|- 1) — E{nu) or E{nul) — E{nu — !)• We have chosen 









the latter definition (Eq. [2]). /iLi is used instead of /XLi+ to control riLi because the former 
is the widely used convention, and because it is harder to control charged particles due to 
the periodic boundary condition adopted. In the literature, Vi is generally averaged over a 
wider range of Li content than in the present work.— 

Ve and Vi are not necessarily equal. A transparent example is a pristine Li metal slab in 
vacuum. Vj=0 by definition. Yet the work function of Li(lOO) is reported to be 2.93 eV.— 
This translates into an electronic voltage of Ve=1.51 V (Eq. [1])! Indeed, if the Li surface 
is covered with a poor SEI that allows e~ transfer but Li+ transfer is slowed to beyond 
experimental timescales - the opposite of what a good SEI does - the lithium slab acts as 
an e~ emitter, the voltage of which must be governed by its E-p- For Li(s) to truly exhibit 
0 V vs. the LIB Li+/Li(s) reference in e~ transfer processes, net surface charge densities 
and electrode-electrolyte interface contributions must be accounted for. Organic solvent 
decomposition can occur on uncharged Li (100) surface under UHV conditions despite the 
relatively high (1.51 V) potential, but in these reactions the metal surface is a reactant and 
Li+ are produced. 

C. Relation between Ve and Vi 

(1) Ve=Vi at equilibrium. Under this condition, Vi controls the Li content. Thus 
Eq. |3] should be interpreted as a self-consistency criterion. /rLi=(*h-1.37 eV-1-1.57 eV) from 
Eqs. [H 131 and the PBE estimate of Epi{s). Identifying $ as /ie, Eq. |2] implies that ^u+ 
becomes a voltage-independent 0.20 eV. This reflects that Li+ is in excess in an infinite 
reservoir. A related point has been made in the fuel cell literature.— In that case 11+ con¬ 
centration at the interface, governed by the pH, is also strongly coupled to the predicted 
voltage. 

(2) When Vi<Ve, and the voltage is held fixed by the potentiostat, and Li+ and e~ sources 
are present, the Li content should spontaneously decrease (as permitted by kinetics) so that 
Vi rises towards Ve- But Vi does not always increase with decreasing nu- Going back to 
the lithium solid example: if the potentiostat is set at 1.51 V in solution, the entire Li slab 
dissolves. If instead an electrically disconnected Li solid is dipped into an electrolyte, some 
Li+ ionizes, leaving a negatively charge, possibly SEI-covered, Li surface. Now it is the 
electronic voltage Ve which decreases to 0 V vs. Li+/Li(s) to coincide with Vi- 
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(3) When Vi>Ve and the voltage is held fixed, the Li content should in general increase 
to lower Vi- This may not occur in all cases. For example, when all available Li+ insertion 
sites are filled in fully lithiated LiCe, the electrode becomes supercapacitor-like. Ve now 
controls the net surface electronic charge,— and V* ceases to matter. Conducting simulations 
under such non-equilibrium conditions are crucial to understanding many phenomena in 
batteries. An example is the kinetics of metal plating^i^i^ at overpotentials - an explicitly 
non-equilibrium phenomenon, a major safety issue in LIB, but also a critical process for 
batteries featuring rechargeable metal anodes. This discussion highlights the primacy of Ve 
when computing voltages. 

Eq. [3] was originally applied to bulk crystalline electrode simulation cells, and V* has 
been identified as the redox potential.— Li should remain in the solid when the applied 
potential is lower than V*. In these special cases of simulation cells without interfaces, the 
Li(s) cohesive energy approach (Vj) is neither consistent nor inconsistent with Vg. This is 
because the electrostatic potential at any point in space is only defined to within a constant 
in the interior of a crystal;— therefore neither Vg nor Ep is well-defined without specifying 
interfaces. Thus metals terminated in different crystal facets exhibit different $’s.— Much 
more significant changes in <F can be induced using adsorbed molecules.— Using DFT-I-U or 
hybrid DFT functional methods, Vi has generally been predicted to be in good agreement 
with measured open circuit voltages in bulk transition metal oxide simulation cells. The 
errors are at most a few hundred meVs.— This suggests that interfacial contributions may 
be small on most battery cathode surfaces. Our work suggests that artificial interfaces can 
be potentially be engineered to give larger, beneficial effects. Molecular adsorption effects 
will be emphasized in Sec. IIVDI 

DFT calculations of interfaces are usually performed in finite-sized, periodically-replicated 
simulation cells kept at overall charge neutrality. In the literature, in DFT calculations 
without ions in the electrolyte (or without any electrolyte at all), the Vi approach has 
therefore assigned a wide range of voltages to LIB interfaces at charge neutrality. In other 
words, a continuous range of potentials-of-zero-charge (PZC) have been assigned by DFT. 
Yet simple “classical” electrodes such as pristine noble metals usually permit one PZC;- at 
voltages away from PZC, the surface is charged. If liquid electrolyte is present, the electrode 
surface can exhibit non-zero charges, compensated by counterions in the electrolyte. In this 
work, we illustrate this effect by adding PF^j" anions under UHV conditions. 
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III. METHODS 


DFT calculations are conducted using the Vienna Atomic Simulation Package (VASP) 
version 5.3^^— and the PBE functional.— Modeling spinel Li 2 ;Mn 204 requires spin-polarized 
DFT and the DFT-I-U augmented treatment^ of Mn 3d orbitals. The U and J values depend 
on the orbital projection scheme and DFT-I-U implementation details; here U — J =4.85 eV 
is chosen to taken from the literature.— The charge state on Mn ions are determined from the 
approximate local spin polarzation Sz reported by the VASP code; |sz|=3/2, 4/2, and 5/2 
are assigned to Mn(IV), Mn(III), and Mn(II), respectively. For Li 3 P 04 and LUCOs, most 
calculations exclude spin-polarization. A few spin-polarized calculations are performed to 
conhrm that this is adequate. A 400 eV planewave energy cutoff and a 3x3x1 Monkhorst- 
Pack grid are applied in all cases. Increasing the cutoff or the /c-grid changes the predicted 
absolute energies by at most 0.1 eV; energy differences, which are the relevant quantities, 
are much smaller. 

Our model systems are asymmetric slabs: Li 2 C 03 basal plane, Li 3 P 04 (010), and 
LiMn 204 (111) slabs on one side of 4-layer thick Au(lll) (Fig. [2]). The standard dipole cor¬ 
rection is applied to negate image interactions in the periodically replicated, charge-neutral 
simulation cells.— An is a metallic conductor and its Up is well dehned. Work functions 
($) are the differences between Up and vacuum levels. Two vacuum-surface interfaces exist 
per asymetric slab, and two vacuum levels, for bare and coated Au(lll) respectively, are 
obtained.— Four layers of An atoms do not give a completely converged <F for uncoated 
Au(lll). Therefore we have shifted the two computed <F of each system by the difference 
between the predicted Au(lll) <h, and the fully converged Au(lll) <F=5.15 eV.— The shift 
is at most 0.1 eV. We also report the net surface dipole density (5), which is the total dipole 
moment of the simulation cell in the direction perpendicular to the surface divided by the 
lateral surface area. The dipole moment is computed as the sum of all charges (including 
electron density dehned on a grid) multipled by their displacement from a user-specihed cen¬ 
ter of the cell, and is reported by the VASP code. 5 should be independent of the cell-center 
in a charge-neutral simulation cell as long as the choice locates the artihcial dipole layer- 
inside the vacuum region. More details on the thin him models, and rationale for crystal 
facets chosen, are provided in the S.I. 
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AA(Li) 

0 

-1 

-2 

-3 

-4 

-1“ 

-U 

Ve 

2.86 

3.26 

3.59 

3.92 

4.19 

3.76 

4.59 

U 

3.07 

3.23 

3.58 

3.65 

NA 

(3.69) 

(4.10) 


TABLE 1: Electronic (Ve) and ionic (Vi) voltages of a Li 3 P 04 ( 010 )/Au(lll) interface cell as the 
number of Li vacancies at the the inner surface of the Li 3 P 04 slab varies. Exceptions: cases a and 
b involve removing a Li from a middle layer or from the outer surface of Li 3 P 04 . Vj is defined for 
incremental Li subtraction, e.g., the “-1” value refers to the difference between the total energies 
of “A” and “-2.” 

IV. RESULTS 

A. Au(lll)/Li 3 P 04 /Vacuum (Fig. [3^ for Summary) 

First we examine the Ve of a non-redox active system: a thin / 5 -Li 3 P 04 (010) solid 
electrolyte him on Au(lll). We start with the stochiometric Li 3 P 04 slab with all Li+ 
sites occupied (AV(Li)=0, Fig. [2^). The work function is $=4.23 eV. This translates into 
Ve=2.86 V vs. Li’'"/Li(s), which is reduced from the bare, charge-neutral Au(lll) surface 
value Ve=3.78 V ($=5.15 V) by 0.88 eV. In the absence of An, the Li 3 P 04 slab is symmetric 
and does not exhibit a net dipole moment. The fact that it modihes the An $ is not due 
to charge transfer between An and Li 3 P 04 . The spatially decomposed electronic density of 
state (DOS) (Fig. 0^) shows that the phosphate Kohn-Sham valence band edge, located at 
about —1.5 eV of the phosphate,— is far from the Fermi level residing on An orbitals. The 
Li 3 P 04 conduction band edge is higher than 5 eV and the vacuum level is at 4.23 eV (not 
shown in the panel). 

Thus the stoichiometric slab must have induced image charges which set up surface 
dipoles. Indeed, the net surface dipole density (5) of the entire simulation cell is predicted 
to be hnite; 5=—0.0048 |e|/A, or —0.023 Debye/A^. The negative sign means that the dipole 
points into Li 3 P 04 from An (111). Using 

AVe = 47r5/(|e|), (4) 

where all quantities are in atomic units, this “small” 5 should yield a large —0.86 V shift 
between bare and phosphate-coated Au(lll). This explains the —0.88 V difference deduced 
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FIG. 4: Kohn-Sham orbitals of Li 3 _ 2 ;P 04 /Au decomposed on to atoms at their coordinates per¬ 
pendicular to the interface. AAi(Li) is (a) 0; (b) —1; (c) —3; (d) —1“; (e) —1^. The vacancies are 
in the phosphate layer closest to Au(lll) (Fig. [5^-b) except for panel (d) (case a, middle of slab, 
Fig. [5b) and (e) (case b, vacancy in outermost layer, not depicted in Fig. [5]). Ep is shifted to zero 
in each case; vacuum levels are not aligned, and are at 4.23, 4.63, 5.29, 5.16, and 6.19 eV in the 
five panels. Au and thin-film-based orbitals are depicted in black and red, respectively. pc=0.007 
is used for determining orbital locations.— 

from work function differences discussed above. DFT calculations give an aggregate 6, that 
includes screening, induced-dipoles, and depolarization effects, not individnal contribntions. 
Hence no effective dielectric constant is present in Eq. 01— 

Under strictly UHV conditions, no net charge can be introdnced. Net bnlk or snrface 
charges lead to inhnite repnlsive energies dne to long-ranged coulombic repulsion. In¬ 
stead, we raise the voltage by removing matched Li^/e“ pairs, as is normally done in DFT 
calculations.—*^ The conhguration space associated with removing mnltiple Li is large. We 
assnme that Li vacancies are weakly interacting, and their formation energetics only depend 
on their positions in the z direction, perpendicular to the surface. Vacancies on the same 
layer of Li 3 P 04 are placed as far apart from each other as possible. With AV(Li)=—4, two 
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(c) 

FIG. 5: Li 3 _a;P 04 /Au: AiV(Li)= (a): -1; (b): -4; (c): -1“. Panel (c) (case a) is associated with a 
Li vacancy in the middle of the phosphate slab. See Fig. [2] caption for color key. 


different configurations with all Li-vacancies right at the interface are indeed found to yield 
total energies within 30 meV of each other. 

Removing one Li (i.e., AA^(Li)=—1) from the Li 3 P 04 layer in contact with An is most 
energetically favorable (Fig. [5^.). It also gives the lowest Ve at this nu content. The DOS 
there reveals no charge transfer between An and Li 3 P 04 (Fig. |Dd). Some spurious Li 3 P 04 
“orbitals” in contact with Au(lll) are above —1.0 eV, but they likely arise from the arbitary 
spatial decomposition scheme; the valence band edge should be taken as —1 eV where at 
least two layers of Li 3 P 04 atoms exhibit the same orbital level.— A net negative charge 
must remain at the Li-vacancy site, compensated by a positive charge on Au(lll) (Fig. [1^). 
This charge separation modifies 5 to —0.0029 |e|/A which should yield a voltage shift of 
—0.53 V from the bare Au(lll) value. Using explicit <F calculations, Ve is found to rise to 
3.28 V, shifted from the Au(lll) Vg by —0.52 V (Table[T]), which is quantatitively explained 
by the change in 5. The phosphate valence band edge is raised by about 0.45 eV relative 
with Ey on the An surface (Fig. |4 ]d). Removing more Li at the phosphate-gold interface 
augments 5 in the positive direction, increases Ve (Table [I]), and raises the valence band edge 
proportionately (Fig. 01:;). Finally, removing a single Li inside the phosphate slab (Fig. 011) 
or at the phosphate-vacuum interface (Fig. 0^) yields a larger Ve than if the vacancy is right 
at the Au(lll) surface (Tabled]). In the latter case, the phosphate valence band edge now 
coincides with Up. 

So far we have not discussed the ionic voltage (Vj) governing the Li content. Table 01 
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shows that Ve exhibits a signihcantly larger slope as a function of AA^(Li) than Vj as long 
as Li are removed from the interior interface. The two voltage dehnitions should intersect 
at only one point. Ve approaches Vi between AA(Li)=—1 to —2. Away from that solid- 
state PZC, equilibrium (Ve=Vi) can be achieved with a net surface charge compensated 
by liquid electrolytes (Fig. [1^, Sec. lIVDj) . just like on noble metal electrodes. Models with 
non-equilibrium interfacial Li content are “wrong” in the sense that they only have transient 
existence. 

Removing a Li sitting deeper inside the phosphate layer (Fig. [T][, [5b) shows signs of 
yielding self-consistency: Ve=3.78 V, Vi=3.69 V (Tabled]). However, the relative free energy 
(AG) of this alternatve AA(Li)=—1 conhguration, evaluated using 

AG = G(AA(Li))-G(0) - AA(Li)/iLi, (5) 

is less favorable than AA(Li)=—3 where all Li vacancies are at the interface (Table dl 
hfth column), at a similar Li chemical potential (/iLi=Vj). In general, if Vj increases upon 
successive Li removal right at the interface, Eq. [5] favors removing multiple Li there, rather 
than removing a single Li further inside the slab - even if the latter yields the same Vi. 
Thus Fig. df should not be a candidate for equilibrium interfaces,— at least under these 
UHV conditions. If a liquid electrolyte is present, it cannot be ruled out that Li+ vacancy 
near the outer surface may be stabilized; we speculate that this could even yield a multitude 
of equilibrium Li”^ conhgurations in Li-containing Elms, and therefore more than one PZC 
at different voltages. 

B. Electrostatics-based Passivation Strategy 

These results emphasize that the interfacial dipole density (5) is one atomic lengthscale 
determinant of applied voltages. From Eq. d] a thin dipole sheet is expected to give rise 
to a step-like voltage drop at the interior interface of Li 3 P 04 (Fig. dfe), which is therefore 
not an inert, passive coating (Fig. db)- At long enough times such that Li’*' can insert and 
relax their conhgurations, this solid electrolyte should absorb much of the voltage drop/rise 
at the sharp metal/electrolyte interface. If a liquid electrolyte exists in the vacuum region 
of Fig. (5] it should experience a reduction in the extreme potential exerted by the cathode 
(Fig. dfe), in effect widening the liquid stability window. Solid Elms coated on anodes should 


15 




AA(Li) 

0 

-1 

-2 

-3 

-4 

-6 

Ve 

3.56 

3.87 

4.03 

4.03 

4.05 

4.41 


3.53 

3.58 

3.26 

3.55 

NA 

NA 


TABLE 2: Electronic (Ve) and ionic (Vj) voltages for a Li 2 C 03 basal plae/Au(lll) interface as 
the number of Li vacancies at the inner surface of the Li 2 C 03 slab varies. No adsorbed molecule 
exists. 

work analogously; in that case, the him, or at least its interface with the anode, must be 
able to incorporate extra Li.— 

The above discussions suggest a new design principle for “smart” artihcial SEI passi¬ 
vating lithium ion battery anodes and high-voltage cathodes. Traditional SEI relies on 
Li’''-conducting but e~ insulating Elms to slow down e~ tunneling from electrodes to organic- 
solvent based electrolytes. They are not designed to affect the voltage experienced by the 
liquid electrolyte, which thermodynamically speaking should still decompose at battery op¬ 
erating conditions. Our proposed scheme focuses on local thermodynamics, not e~ tunneling 
kinetics. 

A fast Li’*' transport rate is crucial for successful implementation of passivation based 
on Li'''-motion induced interfacial dipoles. Another route, not involving Li motion, is to 
incorporate permanent dipoles that can re-orient, such as ferroelectric Elms (Fig. [Tb). Fer¬ 
roelectric surfaces have already been used to manipulate photochemical reactivities.— Most 
ferrorelectrics are not Li’*' conductors. However, electrostatic efiects are long-ranged, and it 
is likely that non-continuous ferrorelectric Elms, lined with gaps that Li’*' can pass through, 
can serve to reduce extreme voltages experienced by liquid electrolytes. 

C. Au/Li 2 C 03 /vacuum (Fig. [3b) for Summary 

The band alignment between Au(lll) and the stoichiometric Li 2 C 03 basal plane (Fig. [2b) 
is similar to that predicted for Au/Li 3 P 04 . In the S.I., the Au(lll)/Li 2 C 03 ( 001 ) interface is 
also shown to be qualitatively similar. The broad correspondence suggests that the predicted 
behavior is universal for non-redox-active Elms. 

The valance and conduction band edges of the AA(Li)=0 basal Li 2 C 03 slab bracket the 
An Fermi level (Fig. [6b), indicating no e~ transfer. V* associated with this Li 2 C 03 slab is 
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predicted to be 3.53 V, very close to Ve=3.58 V there. Thus the PZC somewhat fortuitously 
occurs at stoichiometry. Next, we raise the voltage and try to oxidize Li 2 C 03 .-‘^ by removing 
successive Li’^/e“ pairs. The most favorable Li vacancies are at the Au/carbonate interface 
(Fig. [T^-b). As Li are further removed (Fig. [7J::, Fig. [6j3-c), Ve increases (Table [2]). The 
valence band also rises until the occupied states at the Au/carbonate interface reaches E-p 
(Fig. |6]i). At this point, e~ transfers from LLCOa to Au(lll). Oxidation under UHV 
conditions thus initiate from the Au/Li 2 C 03 interface. The redox reaction pins Ve, which 
ceases to increase with another Li removal (Table E]). This is reminiscent of “Fermi level 
pinning” at metal-semiconductor interfaces by surface/defect states,— except that in LIB 
the “defects” are mobile Li^ vacancies, not static impurities, and their density can be much 
higher. When 6 Li vacancies are created at the interface, Ve increases again. This may due 
to the complete removal of e~ from some surface states. Unlike Li 3 P 04 , the interfacial layer 
of Li 2 C 03 starts to lose structural integrity with removal of multiple Li+ (Fig. [7J::). 

Removing one Li from the middle of the Li 2 C 03 slab (the hole remains on the An sur¬ 
face) is less favorable by 0.25 eV compared to removing Li"^ at the Au/carbonate interface. 
Removing a Li’*' on the vacuum surface (Fig. [Tli) is less favorable by another 0.22 eV, and 
Ve=4.50 V, much higher than the most energetically favorable Li’*' vacancy at the interior 
Au(lll) interface. These conhgurations reveal an inhomogeneity in the energy landscape 
perpendicular to the interface that may strongly affect Li’*' transport^ - especially if Li+ has 
to escape from the interior interface to the outer carbonate surface to the liquid electrolyte. 

D. Au/Li 2 C 03 /EC/vacuum (Fig. [3]c)for summary 

Next we use this Li 2 C 03 example to illustrate the impact of including explicit molecules 
in the theoretical models. A monolayer of adsorbed ethylene carbonate (EC) molecules is 
shown to exhibit profound effect on predicted voltages. Such molecular adsorption might 
be realized under UHV settings. We stress that frozen EC layers at T = 0 K cannot be 
considered adequate approximations of liquids at hnite temperature.— 

Fig. [7^ corresponds to Fig. [Tti except for the adsorbed EC monolayer. The solid him 
atomic structure, energetics of successive Li removal (Vj, Tables |2] and |3]) , and DOS (Fig. |6K 
and Fig. IHk, after shifting Ep to zero), are almost identical with or without EC. This is 
unsurprising because the EC layer neither transfers e~ to nor form a chemical bond with 
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FIG. 6: Kohn-Sham orbitals of Li 2 _a;C 03 /Au decomposed on to atoms at their coordinates perpen¬ 
dicular to the interface, (a); stoichiometric; (b): 2 Li'*' vacancies inside; (c): 3 Li"*" vacancies inside; 
(d): 4 Li'*' vacanies inside. £'f=0.0 eV. The vacuum levels are at 4.93, 5.40, 5.40, and 5.40 eV 
in the four panels. Au and thin-film-based orbitals are depicted in black and red, respectively. 
p^=0.005^ 


FIG. 7: Li 2 _a;C 03 /Au: (a): 1 Li'*' vacancy inside; (b): 2 Li+ vacancies inside; (c): 4 Li'*' vacancies 
inside; (d) 1 Li'*' vacancy outside, (e) 6 EC on stoichiometric Li 2 C 03 /Au; (f) 5 EC, 1 PFg on 
Li 2 C 03 /Au. E is in pink; see Eig. [2] caption for color key. (The hgures are removed due to file size 
limitations - they can be form in the JPCC paper or are available upon requres.) 

Li 2 C 03 . However, this “spectator” EC layer leads to a large Ve decrease. They carry 
sufficiently large dipole moments, even when tilted almost flat in the adsorption geometry, 
to reduce Ve by >1 V compared to vacuum interfaces.— At liquid-solid interfaces, thermal 
fluctuations should modify the 6 value predicted for the monolayer EC in vacuum at T=0 K. 
Nevertheless, liquid EC at hnite temperature has also been shown to change the PZC of the 
LiCe edge plane by more than 1 V.— EC molecules are highly asymmetric; they bind strongly 
to cations, including Li'*' ions at material surfaces,— >21 but weakly to anionic species. The 
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AA(Li),A(PFg-) 

(0,0) 

(-1,0) 

(-2,0) 

(0,1) 

(-1,1) 

Ve 

1.73 

2.07 

2.31 

3.53 

3.45 

Vi 

3.47 

3.61 

NA 

3.57 

NA 


TABLE 3: Electronic (Ve) and ionic (Vj) voltages for the Li 2 C 03 basal plane/Au(lll) interface 
as the number of Li vacancies at the inner surface of the Li 2 C 03 slab varies. The outer surface is 
coated with 6 EC and/or 5 EC and a PEg . 

resulting dipole layer decreases the vacuum Ve on both anodes and cathodes unless the 
surfaces are sufficiently highly charged to force reorientation of the molecules. 

The DOS of Fig. [8 )d may superhcially seem to suggest Li 2 C 03 oxidation can occur at 
a very low Ve=2.31 V. However, Vg predicted with no Li vacancy is already lower than Vi 
computed for removing one Li. The correct interpretation of the hrst 3 columns of Table [3] 
is that the Li content needs to increase in the direction of Ve until all Li-sites are occupied 
(Sec.ini)- It is inconsistent for (Li’'',e“) pairs to be removed from the stoichiometric Li 2 C 03 
slab when Ve=2.31 V. Therefore no oxidation occurs. 

So far we have considered uncharged electrodes. If a liquid electrolyte is present, net sur¬ 
face charges can be compensated by mobile ions in the liquid electric double layer (Fig. [1^). 
Fig. [7f andlHb represent a charged interface after replacing an EC molecule with a PFg .— 
The simulation cell remains charge neutral, and the — |e| charge must be largely compen¬ 
sated by a net positive charge on the Au(lll) surface. This sets up an additional dipole 
moment, which increases <F (Eq. 0]). As a result, Ve rises to 3.53 V, very close to Vi under 
these conditions (Table [3]). The DOS (Fig. [HI:;) shows that, while the PF^ valence band 
remains well below Ep, some Li 2 C 03 orbitals on the outer surface has been raised almost to 
the Fermi level due to the close proximity of the anion. Although the extent of the surface 
charge depletion is difficult to quantify. Fig. [Hb suggests that oxidation may initiate at the 
outer surface of Li 2 C 03 in liquid electrolyte if PFg coordinates to the carbonate surface in 
an inner-shell conhguration (i.e., they are in physical contact), and if the carbonate him is 
thin enough to allow electron transport. 

A thicker carbonate him will increase 5 at the same PFg surface density because of 
the larger charge separation. This means that fewer PF^/ ion near the surface would be 
needed to generate the same voltage increase (Eq. 0]). We speculate that net charges on 


19 










X (A) X (A) X (A) 


FIG. 8: Kohn-Sham orbitals of EC/Li 2 _a;C 03 /Au decomposed on to atoms at their coordinates 
perpendicular to the interface, (a): stoichiometric; (b): 2 Li+vacancies inside; (c): a PF^ replacing 
one EC in panel (a). The vacuum levels are at 3.12, 3.68, and 4.91 eV in the three panels. Au, 
Li 2 C 03 , and EC-/PFfr-based orbitals are depicted in black, red, and blue, respectively. /9c=0.005.— 

thick film surfaces are sparse and can be treated as isolated defects. The anions may also 
be farther away inside the liqnid region and be screened by solvent molecules. As a result, 
the center of negative charge may or may not become more distant from the surface, and 
a quantative estimate of the effect of anion-surface separation on the voltage needs to be 
explicitly computed. Whether thin film-coated electrodes exhibit net charges compensated 
by counterions in the liquid (Fig. [1^), or whether the voltage increase can be accommodated 
entirely by changes in Li^ conhgurations in the solid region (Fig. [!][), depends on which 
option minimizes the Gibbs free energy. 

E. Au/LiMn 204 /vacuum and Au/LiMn 204 /EC/vacuum (Fig. [3]d for Summary) 

In previous work, we have considered EC oxidative decomposition on Li 3 ;Mn 204 (100)^ 
and (111)^ surfaces, in both UHV conditionals'll and in explicit liquid electrolyte.— The 
liquid is found not to strongly affect the rate-determining reaction barrier on the (100) 
surface. The voltage dependence of EC decomposition was not considered there because 
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AA(Li),A(EC),A(PF6) 

(-2,0,0) 

(-3,0,0) 

(-4,0,0) 

(-5,0,0) 

Vg 

3.28 

4.33 

4.71 

5.24 

R 

3.27 

3.58 

4.40 

NA 

AA(Li),A(EC),A(PF6-) 

(-3,1,0) 

(-4,1,0) 

(-5,1,0) 

(-4,1,1) 

Vg 

4.02 

4.23 

4.78 

6.39 

R 

3.86 

4.19 

NA 

NA 


TABLE 4: Electronic (Ve) and ionic (Vj) voltages for the Lia;Mn 204 ( 0001 )/Au(lll) interface as 
the number of Li vacancies at the outer surface varies. 


Li 3 .Mn 204 , like most transition metal oxides, does not have a readily-defined E^- 

In this work, E-p is made unambiguous by putting a thin, reconstructed^ LMO (111) 
slab on a An (111) “current collector” (Fig. [2J::). Removing an e~ from the system, without 
removing an accompanying Li’*', does not immediately lead to a change in any Mn ion 
charge state, strongly suggesting that the e~ is removed from the Fermi level straddling 
An orbitals (not shown). This conhrms that Ep is relevant for voltage estimate via Eq. [H 
Allowing geometry optimization after e~ removal yields a localized hole polaron which is 
now in equilibrium with Ep. 

In contrast to Li 2 C 03 and Li 3 P 04 , successive removal of Li+/e“ pairs from the stoichio¬ 
metric LMO slab occur most favorably from the outer surface (Fig. [9^, see also the S.L), 
not the LMO/Au interface. A Mn(III) turns into a Mn(IV) to accompany each e~ removal. 
Table m indicates that the voltage increases with decreasing surface Li content, and Ve=Vi 
is achieved at AA^(Li)=—2 under UHV conditions. 

Next we add EC to the LMO surface. If the system were a liquid-solid interface, the 
surface Li content, net surface charge, and surface density of EC molecules which directly 
coordinate to available Li ions should be self-consistently deduced using AIMD simulations 
at hnite temperature.— Here we conhne ourselves to T = 0 K and to adding one EC, 
yielding the surface used in Ref. Il7|. AA(Li)=-4 in this model (Fig. |9b), and Ve=4.23 V, in 
reasonably agreement with Vi=4.19 V for this conhguration (Table 0]). The voltage is within 
the experimental LMO operating range. Except for the shift in Vg, and hence in the vacuum 
level, the DOS is almost unchanged from the case where no EC is present iFig. [TOh-bh With 
no anion on the surface, the EDL structure should resemble Fig. [1^. 
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(d) (e) 


FIG. 9: EC and PFg/Lii_2;Mn202/Au. Au atoms are not shown, (a) AA^(Li)=-5; (b) AA(Li)=-4; 
with 1 intact EC; (c) AAi(Li)=-4; with 1 intact EC and 1 PEg ; (d) AA(Li)=-4; with 1 broken 
EC; (e) AA(Li)=-4; with 1 broken EC and 1 PFg . AA(Li) is the number of surface Li vacancies 
compared to the stoichiometric slab. See Fig. [2] caption for color key. 

It is also important to elucidate the EDL structure of a charged, polaron-conducting 
LMO surface, not considered previously.— As in the Li 2 C 03 example, we impose a positive 
charge by optimizing the geometry of a PFg anion on the LMO surface in the overal charge- 
neutral simulation cell (Fig. [91::). The DOS (Fig. flOb ) reveals that PFg retains its excess 
electron and is not oxidized.— A net -l-|e| compensating charge is found in the LMO slab, 
on a second layer Mn ion (“Mn 14”) below the LMO/vacuum interface. This layer contains 
Mn(II) ions, exchanged with Li"*” in the Persson group’s reconstruction,— which may be the 
reason a Mn(III) ion there can lose an e~ most readily. The dipole moment created by the 
Mn(IV)-PF(r charge separation raises Ve dramatically, from 4.25 V to 6.41 V. Comparing 
the DOS with iFig. fTOb) and without iFig. [TObl the PFg , an occupied d-orbital on “Mn 14” 
is indeed seen to be depopulated; its orbital energy level now moves above E-p- The FDL 
in this redox-active material now resembles Fig. [T]i. A lower surface anion density will yield 
smaller shifts Ve but will require a simulation cell with a much larger lateral surface area. 
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F. EC decomposition on LiMn 204 /Au surface 


Finally, we re-examine EC decomposition on the AA(Li)=—4 LMO slab, now with 
Au(lll) present. Fig. [Hll depicts a partially decomposed EC molecule on this LMO sur¬ 
face. A H"*" is transferred to the LMO slab without breaking any C-0 bond in the EC. Its 
energy is -1-0.35 eV above the intact EC conhguration (Fig. [9 ]d), in good agreement with that 
of the reaction intermediate depicted in Fig. 5c of Ref. Il7| where the An slab was absent. A 
subsequent, C-0 bond-breaking step is needed to render EC oxidation favorable.— 

Examining the spin states reveals that a Mn(IV) ion (“Mn 16”) immediately below the 
oxide surface gains an e~ to become a Mn(III). Consistent with this analysis, the DOS shows 
an unoccupied Mn d-orbital (Fig. ITOb for intact EC) becoming occupied (Fig. ITOl i for broken 
EC), moving from above Ep to below it due to polaronic relaxation.— Some orbitals on the 
decomposed EC fragment reside near Ep, suggesting further oxidation can readily occur. 

Regardless of whether an e~ is added to LMO from an oxidized EC, or is removed by 
imposing a -(-|e| charge due to a counter-anion nearby, only Mn redox states change. Charge 
localization on Mn ions, accompanied by polaronic relaxation and a compensating PFg or 
EC+ fragment nearby, ensures that e~ and hole never access the Fermi level which lies on An 
orbitals (Fig. [T0|) . The An current collector only serves to establish Ep, with which polaron 
formation in LMO must remain in equilibrium. Therefore the effect of the applied voltage 
is indirect. Increasing Vi and Ve consistently requires the loss of surface Li’*' (Table llj), 
which leaves the surface with more oxygen ions that are under-coordinated, and increases 
its reactivity towards EC molecules. Adsorbed PFg ions also increases the voltage. 

Since Vi controls the Li content at equilibrium in this voltage range (Li-sites are not 
fully depopulated yet), one can arguably ignore both the An slab and Vg while performing 
simulations of parasitic reaction on such redox-active surfaces.—'^ As discussed above, the 
same reaction intermediate energetics is predicted with or without the An slab. However, 
Ve also controls the PFg surface density. Anions in proximity of the electrode surface 
constitute spatial inhomogeneities or “hot spots” where parasitic or Li’*' insertion reactions 
may preferentially occur. Indeed, when a PF^ is adsorbed near the EC, the deprotonated 
EC intermediate (Fig. [9^) is found to be 0.1 eV more favorable than the intact EC, instead 
of 0.35 eV less favorable when the PFg is absent. Thus volt age-dependent anion adsorption 
can affect interfacial processes. If the cathode material is metallic rather than polaronic. 
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X (A) X (A) X (A) X (A) 

FIG. 10: Kohn-Sham orbitals of EC and PF^/Lii_ 3 ;Mn 204 /Au, decomposed on to atoms at their 
coordinates perpendicnlar to the interface, (a) No adsorbed species, (b) 1 intact EC; (c) 1 intact 
EC and 1 PF^. (d) 1 broken EC. AiV(Li)=-4 in all cases. The vacuum levels are at 6.08, 5.60, 
7.75, and 4.86 eV in the four panels. Au, LMO, and EC-/PFg-based orbitals are depicted in black, 
red, and cyan, respectively. The hlled green and violet circles refer to d-orbitals localized on Mn 14 
and Mn 16 (see text). /9c=0.005.— 


the surface charge is more uniformly distributed,—'^ Ve is more relevant, and qualitatively 
different behavior may be observed. 

V. CONCLUSIONS 

Both electrons and Li+ can move in the surface films covering the electrodes in lithium 
ion battery (LIB), and in the electrodes themselves, in response to voltage variations. This 
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makes LIB rather unique among electrochemical devices. The structures of electric double 
layers (EDL), and the atomic lengthscale manifestations of voltages, become much more 
complex than on pristine noble metal electrodes. 

This work considers Au(lll) slabs coated with Li-conducting solid electrolytes or cathode 
oxide materials under ultra high vacuum (UHV) conditions. The applied voltage is shown 
to be strongly correlated with the surface dipole density at the interface, reminiscent of pho- 
tovoltaics studies. In relatively electrochemically inert thin Elms like Li 3 P 04 and Li 2 C 03 
coated on Au(lll), dipole moments responsible for voltage increase are established by cre¬ 
ating negatively charged Li*^ vacancies at the interface between cathode surfaces and the 
thin solid Elms, compensated with positive charges on the cathode. The phosphate example 
highlights the crucial role of dipole density. The carbonate case emphasizes that organic 
solvent molecules like ethylene carbonate (EC) can dramatically modify the predicted volt¬ 
age even in the absence of charged species (salt). A corollary is that the Poisson-Boltzmann 
theory for liquid electrolytes, which ignores solvent dipole moments, may not quantitatively 
describe battery interfaces. 

With Lia,Mn 204 (LMO) cathode oxide thin Elms coated on Au current collectors, charge 
transfer occurs at the outer LMO surface and involves changes in Mn redox states. The 
Au slab only indirectly determines the snrface Li content and the anion surface density. 
Modeling of electrolyte decomposition on such electrode surfaces arguably does not reqnire 
explicitly including the current collector. However, the voltage-dependent snrface density of 
PFg creates inhomogeneities and hot spots where fast electrolyte decomposition can occnr. 

We also critically examine the deEnitions of “applied voltage” in DFT calculations at 
interfaces. We distingnish two voltages: electronic (Ve), due to electronic motion, and ionic 
(Vi), due to Li+ redistribution. In any atomic conEguration, Ve is well deEned in a metallic 
electrode and should coincide with the experimental voltage imposed via a potentiostat. Vi 
is a self-consistency critierion. At eqnilibrium, the system is pinned by Li+ insertion redox 
reactions, and Ve=Vi. It can be efficiently compnted using the “lithinm metal cohesive 
energy” method widely used in the battery theoretical literature, although most previous 
work neglects the possibility of net surface charges, compensated by ions inside the liqnid 
electrolyte, which may afiect parasitic reactions and Li”^ insertion kinetics. 

But out-of-equilibrium conditions are also critical to many battery-related phenomena 
snch as SEI formation and electroplating. Here Ve, which governs electron content at battery 
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interfaces, is the correct definition of voltage. However, the work function-derived dehnition 
of Ve does not fully account for spatial inhomogeneities. Our solid state modeling consid¬ 
erations form the basis for future AIMD simulations of thin hlm-coated electrode/liquid 
electrolyte interfaces. 

Supporting Information 

A supporting information (S.I.) document provides more details about the model systems 
used, reports calculations on Li 2 C 03 ( 001 )/Au(lll) interfaces and Li 2 C 03 bulk crystal de¬ 
composition thermodynamics considerations, and gives a brief comparison between aqueous 
interfaces and battery electrolyte interfaces. This material is available free of charge via the 
Internet at http: //pubs. acs. org, 
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